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Theoretical studies of the metamagnetism and anomalous phase of SrsRu2 07 have focused on 
the role of van Hove singularities, although much experimental evidence points towards quantum 
criticality having a large effect. We investigate the magnetic and thermodynamic properties of 
systems where magnetic field tunes through such a peak in the electronic density of states. We study 
the generic case of a van Hove singularity in 2D. We see that in combination with the requirement 
of number conservation and interaction effects the peak in the density of states produces several 
interesting phenomena including raising the critical field of the transition above naive estimates, 
altering the relationship between temperature and field scales and creating a distinctive double-peak 
structure in the electronic specific heat. We show that this apparent non- Fermi liquid behaviour 
can be caused at mean- field level by a peak in the density of states. 



There has been much recent interest in the mate- 
rial Sr3Ru207. This displays a metamagnetic transition 
which bifurcates to enclose an anomalous phase where 
the transport properties break the symmetry of the crys- 
tal lattice^. The energetic drive for the metamagnetism 
and anomalous phase formation is generally taken to be 
the presence of a van Hove singularity just below the 
Fermi surface in one of the electronic bands of the ma- 
terial. Such peaks have been shown to reproduce the 
metamagnetic transit iorP^. However, initial studies of 
the material reported a quantum critical endpoiniP in 
the region where the anomalous phase was later discov- 
ered. The role of this critical point in the formation of the 
phase remains unexplored. The region around the phase 
shows signatures which may be attributed to the quan- 
tum critical point, such as a diverging entropy and spe- 
cific heat 8 . In addition the dependence of the critical field 
of the metamagnetic transition on the angle of the ap- 
plied magnetic field, and doping studies 7 have cast doubt 
on the simple picture of a fixed bandstructure where the 
role of the field is simply to Zeeman split the electronic 
spin-species. It is therefore important to identify the dis- 
tinctive signatures of a metamagnetic transition driven 
by a van Hove singularity in order to determine if this is 
indeed the mechanism of the transition and to distinguish 
between quantum fluctuation and bandstructure driven 
effects. 

We will study in detail how the presence of metamag- 
netism and a density of states which varies rapidly near 
the Fermi surface affects the evolution of the system as 
a function of temperature and magnetic field. In partic- 
ular we will focus on the entropy and specific heat of the 
system, quantities which have recently been measured in 
Sr3Ru207 8 . We will study the generic logarithmic di- 
vergence of the density of states due to a saddle-point 
in the dispersion of a 2D system 9 . We will show that 
the peak produces an apparent divergence of C/T and 
unusual behaviour as a function of field. 

We begin by identifying the general mechanisms which 
contribute to the form of the entropy and specific heat 
and predict their behaviour on passing through the meta- 
magnetic transition. We then present the results of a cal- 



culation for the specific heat and entropy from the Stoner 
model for a general density of states. We present the nu- 
merical evaluation of the expressions for the case of the 
logarithmic density of states on a cut through the meta- 
magnetic wing of the phase diagram. These results agree 
well with our general predictions. 



I. FACTORS CONTRIBUTING TO THE 
ENTROPY AND SPECIFIC HEAT 

The shape of the entropy and specific heat curves may 
be deduced from some basic principles. Entropy is a mea- 
sure of the number of available states at the Fermi sur- 
face. At low temperatures it should therefore have the 
same form as the density of states as a function of energy. 
There are several effects which subtly alter this depen- 
dence in the presence of a varying density of states. These 
include spin-splitting due to the external magnetic field, 
interaction induced magnetism, number conservation^ 
and temperature. We will consider each of these in turn. 
We will then examine the specific heat and see that the 
peak in the density of states produces a double peaked 
structure in the specific heat. This is then modified by 
the same effects as the entropy. 

First we consider the entropy, S. S/T is proportional 
to the density of states at the Fermi surface. Upon appli- 
cation of a magnetic field the spin-species' Fermi surfaces 
become split. The entropy is then given by the sum of 
the density of states at two different energies. Since one 
Fermi surface is moved to a lower density of states by the 
spin-splitting this has the effect of compressing the peak 
in entropy around the van Hove singularity as shown in 
panel c) of figure [I] 

We now consider one of the consequences of enforc- 
ing number conservation on a system with a peak in the 
density of states. Applying a magnetic field splits the 
spin-species' Fermi surfaces by an energy proportional to 
the magnetic field. However, by moving a fixed energy 
interval the Fermi surface closest to the peak expands 
to include more electrons than the Fermi surface further 
from the peak loses by contracting (assuming the Fermi 
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surfaces lie below the peak). This means that the overall 
number of electrons has increased. In order to conserve 
number the majority Fermi surface must move towards 
the peak more slowly than the minority Fermi surface re- 
cedes. This results in a slower approach to the peak than 
would naively be expected and therefore a higher critical 
field for the metamagnetic transition. This is shown in 
figure [I] We note that despite this effect the distance 
between the Fermi surfaces remains proportional to the 
field. 

As well as the spin splitting due to the external field 
there is an additional splitting of the Fermi surfaces due 
to the metamagnetism. The effects of this magnetisa- 
tion are sketched in figure [2] When the magnetisation 
is continuous the effect is one of compressing the field 
scale around the crossover so that field tunes through 
the putative entropy curve more rapidly than expected. 
When the transition is discontinuous a range of magneti- 
sation values are 'jumped over' by the transition. This 
removes a slice of the putative entropy curve beginning 
at the critical endpoint and getting wider as tempera- 
ture is decreased and the transition gets stronger. The 
entropy therefore becomes discontinuous at the metam- 
agnetic transition. Since the region removed is around 
the peak in the density of states, it is possible for the 
highest value of S/T to occur at non-zero temperature, 
where the jump is smaller and the Fermi surface sam- 
ples a region closer to the peak. This is similar to the 
susceptibility which is strongly peaked around the criti- 
cal endpoint of the transition, which is expected as this 
point is a second-order phase transition. 

Finally we consider the effect of temperature. This 
broadens the Fermi-Dirac distribution, allowing thermal 
occupation the peak in the density of states across a wider 
range of chemical potential as temperature is increased. 
The peak in entropy therefore becomes broadened as is 
shown in figure ^f). The fact that the Fermi surfaces are 
spin-split means that as temperature is increased each 
will intersect with the peak in the density of states at 
a different temperature, modifying the single peak struc- 
ture of S/T with a broad background from the Fermi sur- 
face furthest from the divergence in the density of states. 

Thermal broadening of the Fermi-Dirac distribution 
also results in the number of electrons changing as the 
density of states is asymmetric about the Fermi level. In 
order to conserve number the chemical potential must 
therefore move as a function of temperature in a similar 
way to that discussed for field. This effect is plotted 
in panel d) of figure [I] where we see it is smaller than 
the effect for field, and that the direction of the motion 
depends on whether the Fermi level is at higher or lower 
energy than the peak in the density of states. 

These effects mean that a comparison of field and tem- 
perature scales is complex. Interactions and number con- 
servation alter the rate at which field tunes into the peak 
in the density of states relative to temperature while the 
seperation of the Fermi surfaces depends not just on the 
magnetic field but on the magnetization. We must there- 



fore take care when associating zero-field features as a 
function of temperature with low temperature features 
as a function of field. These details of these effects will 
be studied in more depth in a later publication. 

The low temperature specific heat of non-interacting 
fermions normally follows the density of states, however 
in the case of a peak in the density of states near to the 
Fermi level and this approximation fails. We will consider 
how temperature, magnetic field, and the inclusion of 
interactions change this result. 

Raising temperature has the effect of broadening the 
Fermi-Dirac distribution, the range of states occupied 
around the chemical potential is therefore increased. 
Starting with the chemical potential on either side of the 
density of states peak and raising the temperature will 
allow thermal occupation of the states under the peak. 
As shown in figure [3] this increase is most rapid on ei- 
ther side of the peak. Specific heat is proportional to the 
temperature derivative of the entropy and will be largest 
in these regions. We therefore expect a structure which 
is logarithmic at zero temperature with a peak which bi- 
furcates and broadens as temperature increases. 

The magnetic field and interactions have the same ef- 
fect as discussed for the entropy - a compression along 
the field axis and the removal of a 'wedge' of field values 
due to the first order transition. The rate of approach 
to the transition will be altered by number conservation 
as discussed. This will lead to C/T deviating from the 
symmetric logarithmic structure as we will see in section 

eh 

We will now consider how to calculate these quantities 
explicitly within the Stoner mean-field theory. 



II. DERIVATION OF EXPRESSIONS FOR 
ENTROPY AND SPECIFIC HEAT 

We calculate the entropy and specific heat from the 
free energy of the Stoner model. The effective chemical 
potential for the spin-species is determined by number 
conservation. Number conservation is enforced by requir- 
ing n a = | + am where a = ±1 labels the spin-species, 
the total number of electrons n = + is constant 
and the magnetization is given by 2m = — n^. The 
effective chemical potential \i G is determined implicitly 
from the number n a , 



de p(e)f(e - 



(1) 



-l . 



where f(e — p a ) = (l + exp (e — p))) " is the Fermi- 
Dirac distribution, we have set &b = 1, and p(e) is the 
density of states. 

The free energy in the Stoner model is: 



p a n 



(2) 
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FIG. 1. a) When splitting the spin-species' Fermi surfaces with a magnetic field, number conservation forces a non-linear 
dependence of the chemical potential on field. The black line is the zero- field filling and the Fermi surfaces are split by 2H. The 
dotted lines show the case where total number is allowed to vary. Solid lines show the situation where number conservation 
has been enforced. The Fermi surfaces are still split by 2H but have both moved to lower energy so that there are the same 
number of electrons as before the splitting, b) Solid lines are chemical potentials for the up- and down-spin electrons and 
their average. We see that the Fermi surfaces move with field at different speeds and are slowed in their approach to the peak 
compared to the non-number conserving case shown by the dashed lines, c) 'Entropy' curves constructed from the sum of the 
density of states at each Fermi surface. Blue (solid line): If field tuned a single Fermi surface through the density of states. In 
this case the entropy would mirror the density of states. Purple (dashed line): Taking into account spin-splitting. The peak is 
compressed due to one spin species sampling a lower density of states. Red (dotted line): Number conservation results in a a 
higher field to tune through the singularity, d) Chemical potential as a function of temperature. Number conservation results 
in \i changing as a function of temperature. 



where g is a constant representing the interaction 
strength and H is the magnetic field. From the require- 
ment that the free energy is a minimum we obtain a self- 
consistent equation for the magnetization and magnetic 
susceptibility 3 , 



H /i t (n, m) - /i;(n, m) - 2gm, 



fde p(e)dj(e - p 



2g. 



(3) 



(4) 



S = J2\fdep(e)ln(l + e-^) 
+ T J de p(e)d T \n (l^e'^ 
C = E Jdeep(e)d T f(e-^) 



(2gm + H) drn^. 



(5) 

The temperature derivatives of the chemical potential are 
given by 



d T Po 



fde 



2gfde E a e -^ 



1 _ M 

± T 
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The entropy is defined by S = — drF\ n H and the specific where 
heat as C = —T d^F\ „. These are evaluated with the 

1 \n,H 

condition that total number is conserved. This condition 
is encoded in the behaviour of the chemical potentials, 
giving a non-trivial form for <9x/V- The evaluation of 
these derivatives is straightforward but lengthy. Here we 
give the results of these calculations: 



P(e) 



(e-Mcr) 

e t 



(6) 



(7) 



The specific heat and entropy may be calculated for 
any n, H and T from equations (|3]),([5]) and These 
expressions give the magnetic transitions of the Stoner 
model, although the location of the first-order transition 
must be determined by minimisation of the free energy 
as it is not uniquely determined by ([3| which becomes 
mult i- valued in the region of the first-order transition. 
We will now evaluate these expressions for magnetization, 
magnetic susceptibility, entropy and specific heat for the 
generic logarithmic density of states. 
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FIG. 2. On the top row we sketch the magnetisation and 
entropy in the absence of the metamagnetic transition. The 
values which will be removed by the transition are shaded. On 
the bottom row the first order transition removes a 'wedge' of 
magnetization values, this removes a section from the putative 
entropy curve. 
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FIG. 3. Temperature broadening of S/T with spin-splitting 
present. Blue, low temperature, red, high temperature. On 
the right the difference between these curves reveals the 
double-peak structure. Note that due to a multiplicative fac- 
tor of T the central region of the peak does not have a negative 
S. 



III. DENSITY OF STATES WITH 
LOGARITHMIC PEAK 



Figure Elgives the results of evaluating (|3| for magne- 
tization, for magnetic susceptibility and ([5| for en- 
tropy and specific heat, for a cut through the metamag- 
netic wing in the H, T plane of the phase diagram with 
the logarithmic density of states ([sj) . These plots are in 
good agreement with the anticipated results. Magnetiza- 
tion has the familiar first-order transition at low temper- 
ature which becomes continuous at a critical endpoint. 
The susceptibility shows a line of maxima tracing the 
magnetic crossover above the critical temperature of the 
transition. There is a peak at the critical endpoint where 
the phase transition is second order. The entropy has a 
temperature-broadened peak reflecting the form of the 
density of states with the position and symmetry of the 
peak shifting due to the effects of the metamagnetism and 
number conservation. Specific heat shows the expected 
double-peak structure but with very asymmetric peaks. 
This is due to the fact that spin-splitting means there 
is a contribution from the minority Fermi surface which 
is different on each side of the transition, also the mag- 
netic transition means that the rate of spin-splitting on 
the high and low magnetization sides is different. Other 
signatures of the metamagnetic transition are observed 
in the field dependence, such as the discontinuity at the 
first-order transition and the maximum of S/T at non- 
zero temperature. 

As would be expected for this form of density of states 
the specific heat has a dominantly logarithmic depen- 
dence on temperature immediately above the metam- 
agnetic transition. This dependence is modified by the 
presence of spin-splitting and number conservation as be- 
comes apparent at higher temperatures than those plot- 
ted in figure [4] and will be discussed in more detail in a 
later publication. 

The ratio between the zero field peak in the specific 
heat and the critical field of the transition can be tuned 
by varying the interaction strength. This illustrates 
the point made about comparing field and temperature 
scales. 



Using the previous results it is possible to evaluate the 
magnetization, entropy and specific heat as a function 
of band filling, magnetic field and temperature, for any 
given density of states and interaction strength. We will 
study a logarithmically divergent density of states, as 
produced by saddle points in the electronic dispersion. 
This model density of states is given by 



P(e) = ^ln 



W 



(8) 



where the bandwidth is 2W . The density of states di- 
verges at e = e c . In the following we will take the inter- 
action strength to be g = 0.31^, the form of the phase 
diagram remains the same if this value is altered. We 
choose to look at a filling which is below the van Hove 
point and use magnetic field to tune the system through 
the metamagnetic transition. 



IV. CONCLUSIONS 

A van Hove singularity is known to cause a metam- 
agnetic transition within the Stoner theory. As well as 
the magnetization and magnetic susceptibility this fea- 
ture in the density of states and the associated magnetic 
transition have distinctive features in the entropy and 
specific heat of the system, resulting in an apparently 
non-Fermi liquid type behaviour. We have calculated 
this behaviour for the generic case of the 2D van Hove 
singularity. S/T reflects the form of the density of states, 
being a temperature broadened peak. The divergence of 
this peak as temperature is lowered is cut off by the meta- 
magnetic transition. C/T shows an asymmetric double- 
peaked structure, also appearing to diverge as a function 
of temperature at the critical field. We have shown that 
the requirement of number conservation alters the rate at 




FIG. 4. From top to bottom the magnetization, magnetic susceptibility, entropy and specific heat as we cross a metamagnetic 
transition in the Stoner model with a logarithmic peak in the density of states, g — 0.3VK, filling fraction is 0.42. The results 
are shown both as 3D plots and gradient plots for clarity, jagged edges in the figures are due to innaccuracies in the plotting 
method. 
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which magnetic field tunes through the peak in the den- 
sity of states meaning that the field dependence of these 
quantities is not simply a constant sweep through the 
density of states. At the first-order transition both en- 
tropy and specific heat jump. This discontinuity vanishes 
at the critical point but the magnitude of the entropy on 
the high-field side of the transition is not monotonic in T, 
having a maximum at non-zero temperature. These re- 
sults suggest that careful consideration should be made 
of the effects of the bandstructure before trying to ex- 
tract information about the quantum critical properties 
of magnetic transitions as even the mean-field analysis 



produces striking features. This study was motivated by 
the material Sr3Ru207 but the results are applicable to 
any system where the Fermi surface lies close to a van 
Hove singularity. 
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